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ABSTRACT 

The goal of our serial papers is to examine the evolutionary connection be- 
tween AGN and star formation in its host galaxy in the partially obscured AGNs 
(i.e., Seyfert 1.8 and 1.9 galaxies). Taking advantage of these galaxies, the prop- 
erties of both components can be studied together by direct measurements. In 
this paper, we focus on the broad-line composite galaxies (composite AGNs) 
which are located between the theoretical and empirical separation lines in the 
[N II] /Ha vs. [O III]/H/3 diagram. These galaxies are searched for from the com- 
posite galaxies provided by the SDSS DR4 MPA/JHU catalogs. After re-analyze 
the spectra, we perform a fine classification for the 85 composite AGNs in terms 
of the BPT diagrams. All the objects located below the three theoretical separa- 
tion lines are associated with a young stellar population (<lGyrs), while either 
a young or old stellar population is identified in the individual multiply-classified 
object. The multiply-classified objects with a very old stellar population are lo- 
cated in the LINER region in the [O I] /Ha vs. [O III]/H/3 diagram. We then 
consider the connection between AGN and star formation to derive the key re- 
sults. The Eddington ratio inferred from the broad Ha emission, the age of the 
stellar population of AGN's host as assessed by D n (4000), and the line ratio 
[O I] /Ha are found to be related with each other. These relations strongly sup- 
port the evolutionary scenario in which AGNs evolve from high L/Lsdd state with 
soft spectrum to low L/L^d state with hard spectrum as young stellar popula- 
tion ages and fades. The significant correlation between the line ratio [O I] /Ha 
and _D n (4000) leads us to suggest that the line ratio could be used to trace the 
age of stellar population in type I AGNs. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: evolution 
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1. Introduction 

1.1. Composite galaxies 

The Baldwin-Phillips- Terlevich (BPT) diagrams (Baldwin et al. 1981) are commonly 
used as a powerful tool to classify narrow emission-line galaxies according to their domi- 
nant energy sources. The diagrams are then refined by Veilleux & Osterbrock (1987). In 
the diagnostic diagram [O III]A5007/H/9 vs. [N II]A6583/Ho;, the star- forming sequence 
shows a curved concentration extending from up- left to bottom-right, while AGNs (includ- 
ing LINERs, Low Ionization Nuclear Emission-line Regions) concentrate on the right part of 
the diagram. A theoretic demarcation curve determining the upper limit of the theoretical 
stellar photoionization model is proposed by Kewley et al. (2001), and an empirical curve 
dividing the pure star-forming galaxies recently by Kauffmann et al. (2003a). The region 
between the two curves is occupied by the so-called composite objects (hereafter compos- 
ites for short) whose spectra are believed to contain significant contributions from both star 
formation and AGN (Seyfert or LINER). 

As a fundamental, tight relationship, the well-established relation between the mass of 
supermassive black hole (SMBH) and velocity dispersion of the bugle in which the SMBH 
resides is served as a clue of co-evolution of AGN and star formation in the bugle (e.g., 
Magorrian et al. 1998; Gebhardt et al. 2000; Granato et al. 2001; Tremaine et al. 2002; 
Ferrarese & Merritt 2001; Greene & Ho 2006). At first, the issue of co-evolution could 
be studied in phenomenology. In fact, there are many observations of both imaging and 
spectroscopy that show the evidence of co-existence of AGN and circumnuclear star formation 
(e.g., Canalizo & Stockton 2001; Cid Fernandes et al. 2001, 2005; Hines et al. 1999; Schinerer 
et al. 1998; Heckman et al. 1997; Gonzalez Delgado 2002; Zuther et al. 2007; Sarzi et al. 
2007; Asari et al. 2007; Davies et al. 2007). 

We should next answer the questions: What is the physics that links star formation and 
AGN phenomena? How AGN co-evolves with star formation in details? (e.g., Ho 2005a). So 
far, there are two distinct scenarios that connect AGN and star formation phenomena: 1) 
starburst occurs first and triggers AGN's activity latter (e.g., Weedman 1983); 2) in contrast, 
AGN's activity occurs before the beginning of starburst (Goncalves et al. 1999 and references 
therein). If the hypothesis of co-evolution of AGN and the bulge of its host is correct, then 
whatever AGN occurs before/after starburst, the composites might represent an important 
state in galaxy evolution, and are valuable for examining the issue of co-evolution in details. 
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1.2. A new approach by Seyfert 1.8 and 1.9 galaxies 

In order to study the co-evolution of AGN and star formation, properties of both AGN 
and its host galaxy are required to be examined carefully and simultaneously. However, 
discrimination between type I AGN and its host galaxy is a difficult task, because both 
the AGN's continuum and broad lines mask the stellar features in optical wavelength. The 
problem of AGN-host separation could be partially solved by several techniques listed as 
follows: 

• Off-nuclear spectroscopy: By detecting the Balmer absorption lines and emission lines 
of H II region in the surrounding of QSOs, Canalizo & Stockton (2001) examined 
properties of circumnuclear stellar populations in the hosts of nine local QSOs. The 
estimated ages range from a few Myrs old to post-starburst. In a similar way, Nolan 
et al. (2001) found that the off-nuclear (~5") stellar populations of optically selected 
QSOs are generally dominated by old ones (~8-14Gyrs). Relatively old stellar pop- 
ulations (~0.1-2Gyrs) are also frequently detected in the off-nuclear regions in a few 
powerful radio-loud AGNs (Tadhunter et al. 2005). 

• AGN+post starburst: In this case the star light is easily discriminated in the 1-D spec- 
troscopy because of the prominent Balmer absorptions of A-type stars. Brotherton et 
al. (1999) identified a 400Myrs old instantaneous starburst in quasar UN J1025-0040. A 

~ lOOMyrs old post-starburst is identified in NLSl-like type I AGN SDSS J085338.27+033246.1 
(Wang & Wei 2006). Wang et al. (2004) identified an evident Balmer jump in type I 
AGN SDSS J022119. 84+005628.4. 

• Spectroscopy of obscured AGNs: By analyzing the spectra of narrow-line AGNs from 
Sloan Digital Sky Survey (SDSS; York et al. 2000), Heckman et al. (2004 and references 
therein) and Kewley et al. (2006) found that most of the accretion-driven growth of 
mass of SMBH occurs in the galaxy associated with a relatively young stellar population 
(see also in Heckman & Kauffmann 2006). 

• Infrared color: By analyzing a sample of IRAS-selected type 1.5 AGNs, Wang et al. 
(2006) found that the essential Eigenvector 1 (El) space (Boroson & Green 1992, 
hereafter BG92) is related with the age of stellar population of the host as assessed by 
infrared color a(60,25). 

• Deconvolution of 2-D spectra: By performing the MCS deconvolution method on 2-D 
spectra, Letawe et al. (2007, and references therein) recently separated the individual 
spectra of the quasar and of the underlying host galaxy in a sample of 20 low-redshift 
luminous quasars. Quite young stellar populations (~ 1 — 2 Gyrs) are subsequently 
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identified in a sample of 18 quasars by a semi-analytic decomposition (Jahnke et al. 
2007). 

Although the understanding of the co-evolution of AGN and star formation has been 
remarkably developed over the past few years by using these approaches, their intrinsic dis- 
advantages prevent the systematical and extensive investigation on the issue of co-evolution. 
We refer the readers to Letawe et al. (2007) for a comment of the drawbacks of these 
approaches. In the present study, taking advantage of Seyfert 1.8 and 1.9 galaxies (i.e., 
partially obscured AGNs), a new approach is adopted by us to directly obtain properties of 
both central AGN and circumnuclear stellar population systematically and simultaneously. 
Seyfert 1.8 galaxy is defined as one whose spectrum shows strong narrow components, and 
very weak but still visible broad components of Ha and H(3. No broad components of 
Balmer lines except a weak broad component at Ha can be identified in Seyfert 1.9 galaxy 
(Osterbrock 1989). In these two types of AGN, due to the orientation effect (e.g., Antonucci 
1993), their optical continuum emitted from the central accretion disks is blocked by a 
torus, and the star light from the bugle is much easier to be detected. In addition, the broad 
emission lines shown in the spectrum, rather than some tracers, allow us to derive AGN's 
properties directly. 

As the first step, by deriving the properties of both AGN and stellar population from 
the optical spectra taken by SDSS, we focus on the nature of the composites with broad 
Balmer emission lines (hereafter composite AGNs for abbreviation), and then study their 
implications on the issue of co-evolution in this paper. The issues related with other types 
of galaxies (e.g., Seyfert and LINER) will be presented in subsequent paper(s). At present, 
SDSS provides a large and uniform spectroscopic sample of galaxies to make such study 
available. The paper is organized as follows. The sample selection and data reduction 
are described in Section 2. The physical properties of both AGN and stellar population 
are derived in Section 3. Basing upon the BPT diagrams, a further classification of the 
composite AGNs is shown in Section 4. We investigate the AGN-host connection in the 
composite AGNs in Section 5. The A cold dark matter (ACDM) cosmology with parameter 
h = 0.7, Qm = 0.3, and = 0.7 (Spergel et al. 2003) is adopted throughout the paper. 

2. Sample Selection and Data Reduction 

2.1. The Sloan Digital Sky Survey 

The SDSS is designed to eventually provide an imaging and spectroscopic survey of 
one-quarter of the entire sky (e.g., York et al. 2000). The survey is carried out by a 
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dedicated 2.5-meter wide-field (3 ° ) telescope at Apache Point Observatory. The telescope 
is equipped with two fiber-fed spectrographes and a mosaic CCD camera. Each spectrum 
is taken with a 3" diameter fiber aperture, and covers a wavelength range from 3800A to 
9200A at the observe-frame. The spectral resolution varies from 1850 to 2200, corresponding 
to (Tint ~ 65km s _1 . The spectrophotometric calibrations with accuracy of about 20% are 
performed by observing subdwarf F stars in each 3 ° field of view. Tow pipelines, spectro2d 
and specBS are developed to reduce the raw spectra (Glazebrook et al. 1998; Bromley et al. 
1998). 



2.2. Selection of Seyfert 1.8 and 1.9 galaxies 

We select a sample of the composite AGNs from the SDSS-Data Release 4 (Adelman- 
McCarthy et al. 2006) MPA/JHU catalogs (Heckman et al. 2004; Kauffmann et al. 
2003a,b,c; Tremonti et al., 2004)0. The catalogs provide spectral and physical properties 
for 567,486 individual narrow emission-line galaxies, including Seyferts, LINERs and star- 
forming galaxies. 

First of all, in order to ensure that our results are not distorted by the spectra with poor 
quality, we only consider the galaxies whose spectra have median signal-to-noise ratio per 
pixel of the whole spectrum S/N > 20. Secondly, the star- forming galaxies which are under 
the demarcation line given in Kauffmann et al. (2003a) are excluded by simply adopting the 
emission line intensities provided by the MPA/JHU catalogs. The extinction is not corrected 
here since both the line ratios [O III]/H/3 and [N II] /Ha are not sensitive to dust extinction. 
There are total 27,341 galaxies that fulfill these criterion. An object is then classified as a 
composite if 

log([OIII]/H/3) < 0.61/(log([NII]/Ha) - 0.05) + 1.3 (1) 

and 

log([OIII]/H/3) > 0.61/(log([NII]/Ha) - 0.47) + 1.19 (2) 

Altogether, there are 12,343 objects that could be classified as the composites with median 
S/N ratio lager than 20. 

An automatic routine is then developed by us to reduce the raw data and to select 
objects with broad emission lines attributed to central AGN. At the beginning, the Galactic 
extinction is corrected for each spectrum by assuming an Ry = 3.1 extinction law (Cardelli et 
al. 1989). The color excess E(B—V) for each spectrum is taken from the Schlegel, Finkbeiner 



These catalogs can be downloaded from 



http: / /www. mpa-garching.mpg.de/SDSS / 
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and Davis Galactic reddening maps (Schlegel et al. 1998). Each of the extinction-corrected 
spectra is then transformed to the rest frame, along with the /c-correction, according to the 
corresponding redshift provided by the SDSS pipelines. 

The contribution of star light is then separated from each spectrum by the principle 
component analysis (PCA) method (e.g., Yip et al. 2004; Li et al, 2005; Hao et al. 2005). A 
library of stellar absorption spectra is built by applying the PCA technique to the standard 
single stellar population (SSP) models developed by Bruszual & Chariot (2003, hereafter 
BC03). The models have a resolution of 3 A. The optical spectral evolution is provided for 
six different metallicities. At each metallicity, there are total 221 spectra at the unequally 
spaced time steps from the age of O.lMyrs to 20 Gyrs. The first seven eigenspectra are used 
to model the starlight component simply through their linear combination. The component 
of AGN's continuum is ignored in the separation in most of the cases, since the contribution 
of the observed continuum of AGN is negligible in their spectra (e.g., Kauffmann et al. 
2003a), except in 7 objects. In summary, in order to appropriately separate the starlight, 
our template contains the seven eigenspectra and a Galactic extinction curve (Cardelli et 
al. 1989) for most of the cases. Because only the instrumental resolution is taken into 
account in the models of BC03, the template is convolved with a Gaussian function to 
match the velocity dispersion in each object. The velocity dispersion is determined by the 
cross-correlation method before the fitting and subtraction. For the cases that show emission 
of both broad H/5 line and Fe II complex of AGN, both an additional power-law continuum 
with form F u oc u~ - 5 (Richards et al. 2006; Elvis et al. 1994) and a template of the Fe II 
complex are required to model the spectra. We simply adopt the empirical template of the 
Fe II complex given in BG92. The line width of the template is forced to be identical to that 
of H(3 broad component (e.g., BG92). 

A x 2 minimization is used to model the starlight component in each spectrum. The 
fitting is performed over the rest frame wavelength range from 3700A to 7000A, except 
for the regions around the strong emission lines: e.g., Balmer lines, [O III]AA4959, 5007, 
[N II]AA6548,6583, [S II]AA6716, 6731, [O II]A3727. The possibly broad components of 
Balmer lines are also considered to determine the emission-line free regions. As an illus- 
tration, the modeling of the starlight component is shown in the left column of Figure 1 for 
three typical cases. In each panel, the upper and lower curves shows the starlight-subtracted 
(emission-line) spectrum and modeled starlight spectrum, respectively. 

The emission-line spectra are subsequently used to search for the composite AGNs, 
according to the existence of a high velocity wing of broad Ha component at its blue side. 
The red wing is not used in the sample selection because of the superposition of the strong 
[N II]A6583 emission line. The existence of a broad wing is accepted if the flux ratio F w /a c > 
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3, where F w is the flux of line wing averaged over the wavelength range from 6500A to 
6530A at the rest frame, a c the standard deviation of flux within the emission-line free 
region between A5980 and A6020. Iterative tests on a typical small sample taken from 
SDSS are performed to determine the wavelength ranges quoted above. The automatically 
selected composite AGNs are then inspected by eyes. Finally, our sample totally contains 
85 composite AGNs with broad Ha emission line. 



2.3. Emission line measurements 

The emission lines are measured once the starlight component is subtracted. We model 
and measure emission lines of the 85 composite AGNs by the SPECFIT task (Kriss 1994) in 
the IRAF packag^E In each spectrum, each emission line is fitted by a set of Gaussian com- 
ponents. The intensity ratios of the [O III] and [N II] doublets are fixed to their theoretical 
values (e.g., Dimitrijevic et al. 2007). In Figure 1, the emission-line modeling in H/? region 
is schemed in the middle column, and the modeling in Ha region in the right column. Direct 
integration is used to measure the flux of [O I]A6300 emission line in each spectrum through 
the SPLOT task. The uncertainties of the emission-line measurements are dominated by the 
subtraction of the starlight component, especially for the faint H/3 emission. Similar as Asari 
et al. (2007), a wide (~ 200A) absorption trough in the continuum around H/3 is frequently 
identified in our starlight-subtracted spectra, especially in the galaxies associated with an 
old stellar population, which leads to a 2%-7% underestimate in H/? flux if the continuum is 
fixed to be zero. This small mismatch is probably caused by the calibration in the STELIB 
library (Asari et al. 2007). We ignore this small underestimate in this paper because rel- 
atively young stellar populations are frequently identified in our sample (see Figure 5 for 
details). 

The results of the emission-line measurements are listed in Table 1. Column (1) lists 
the identification of each object listed in our sample, and Column (2) the corresponding 
redshift. For each object, the line ratios of [N II]/Ha, [S II]AA6726, 6731/Ha, [O I]/Ha 
and [O III]/H/3 are listed from Column (3) to Column (6). Column (9) and (12) shows 
the intrinsic luminosity of Ha broad component (Lh«) and that of [O III]A5007 (L[oiiij), 
respectively. Both luminosities are corrected for local extinction. The extinction is inferred 
from the narrow-line ratio Ha/H/3 for each object, assuming the Balmer decrement for 



2 IRAF is distributed by the National Optical Astronomical Observatories, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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standard case B recombination and Galactic extinction curve with Ry = 3.1. 



3. Deriving Physical Properties 



In this section we derive physical properties of both AGN and stellar population basing 
upon the reduction and measurements described above. 



There are compelling evidence suggesting that both Eddington ratio (L/Le) and virial 
mass of SMBH (Mbh) are basic parameters determining the properties of individual AGN 
(e.g., Boroson 2002; BG92; Sulentic et al. 2000; Grupe 2004). In this study, we attempt to 
directly obtain the physical properties of individual central AGN according to its broad emis- 
sion lines. Basing upon the calibration of the measurements of the reverberation mapping 
methods, the empirical relationships that are commonly used to estimate L/L E and Mbh in 
typical type I AGNs have been well established in recent years (e.g., Kaspi et al., 2000, 2005; 
Peterson et al. 2004; Vestergaard et al. 2002; Wu et al. 2004). Both velocity of broad H/3 
emission line and luminosity at 5100A (or luminosity of broad H/5 emission) emitted from 
the central AGN are required in these calibrations, which are however not feasible in our 
sample both because of the dilution of AGN's luminosity by the host's starlight and because 
of the non-detectable broad H/3 emission in a majority of cases. Greene & Ho (2005) recently 
extended the estimations of both bolometric luminosity and M B h to broad Ha line 



The Ha emission line provides much more robust estimations of both L/L E and Mbh in the 
cases where the contamination caused by host and/or non-thermal jet radiation is present. 
The relationships allow us to estimate both L/L E and M B h in the composite AGNs. For 
each object, the estimated values of M B h and L/L E is listed in Column (10) and (11) in 
Table 1, respectively. 

According to the unified model of AGN (Antonucci 1993), the existence of the broad 
emission lines in AGN's spectrum depends on the orientation of the axis of tours with 
respect to the line-of-sight of observer. In order to assess how much the orientation effect 
affects Lu a (and hence L/L E dd and Mbh) that is estimated from the broad Ha line, the 



3.1. Deriving properties of AGN 




(3) 
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£[oni] vs - Lu a diagram is shown in the Panel D in Figure 7. A strong correlation is found 
between L[om] and L Ha . [O III] emission is isotropic in AGN because the line is emitted 
only from the optically thin region that is at large enough distance from the central black 
hole. The isotropy of [O III] emission is also supported by the correlation between [O III] and 
orientation-independent [O II] emission (Kuraszkiewicz et al. 2000), although the isotropy 
has been questioned in some radio-loud AGNs (Baker 1997; di Serego Alighirei et al. 1997). 
The tight correlation shown in the Panel D in Figure 7 therefore indicates that the variation 
of L Ha in our sample is not entirely due to the orientation effect, and that L Ha is a good 
measure of AGN's power in our sample. 

A simple least-square linear fit within a wide range of luminosity gives an empirical 
relationship between Lpm] and L Ha 

log ( t mi] t ) = (1-518 ± 0.06) + (1.078 ± 0.07) log ( —J^ T ) (5) 

& VlO 40 erg s- 1 ; V ; V ; & V1° er g s / 

Incorporating with Eq. 3 yields a relationship between continuum luminosity L 5W0 and 
Ljoiiij: L 5100 /10 42 erg s _1 1.5(L[oni]/10 40 erg s -1 ) 8 . Taking into account of the dynamical 
range of L[oiii] from ~ 10 40 to ~ 10 42 erg s _1 , an average ratio is -^sioo/^oni] ~ 100. 
Although the ratio is three times lower than the mean value derived in Kauffmann et al. 
(2003a) and in Heckman et al. (2004), it is reasonable because of the large scatter in the 
correlation between Lpm] and L 5WQ . There are two possible explanations for the discrepancy 
of the L 5W0 /L[om] ratio: 1) The luminosity of the broad Ha line (and hence L/L-^dd and Mbh) 
might be underestimated due to the obscuration caused by torus in the frame of the unified 
model of AGN (Antonucci 1993). It is noted that this systematic shift would not seriously 
affect our final statistic trends, both because of the tight correlation between L Hq and Ljoiii] 
and because of the weak dependencies of L/Lgcid and Mbh on Lh<* (i-e., L/L-^dd ^ an( l 
Mbh oc L^ 5 , see Eqs. 3 and 4). 2) In additional to the AGN's contribution, star formation 
can also contributes to the [O III] luminosity. In fact, our sample is mainly composed of the 
composites in which the narrow emission lines are likely to be mainly ionized by hot stars 
(e.g., Kewley et al. 2006, and see below). 

We estimate the upper limits of the underestimates of both L/L-^dd and Mbh as follows. 
Assuming the underestimate of Lh« is entirely caused by the obscuration of torus and by the 
reddening in BLR, the observed broad Ha luminosity L' Ha could be written as L' Ha = ftfeLna, 
where L Hq is the intrinsic Ha luminosity, f t and f e represents the decreases of the line 
luminosity due to the obscuration and reddening, respectively. If the three times lower in 
the L 5 ioo/£[oiii] ratio is caused by the obscuration, f t is inferred to be ~ 0.28 according 
to the relation L 5 i o oc L^. The factor f e could be estimated from the 19 objects for 
which the broad H/9 is present. Recent studies found the mean value of the broad-line 
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Balmer decrements (Ha/H/3) in typical type I AGNs is significantly close to the value of 3.1 
predicted by the case B recombination (e.g., Greene & Ho 2005; Vanden Berk et al. 2001; 
Dong et al. 2007). However, large Balmer decrements (7-10) are observed in Seyfert 1.9-like 
galaxies (e.g., Osterbrock 1981). Given the Balmer decrement in the case B recombination, 
only upper limit of f e could be obtained in our sample, both because of the orientation effect 
and because of the fact that higher density and/or higher ionization parameter is required to 
produce H/3 than Ha. The mean value of the broad-line Balmer decrement (Hq;/H / 9)b is 5.97 
in current sample. Excluding the reddening inferred from the narrow-line Balmer decrement 
(Ha/H/3) N = 4.13, we obtain f e ~ 0.47. Finally, the upper limit of the underestimate of 
L/L Edd is estimated to be ~ 50%, and that of M BH to be ~ 70%. 



3.2. Deriving age of stellar population of host 

Both 4000A break (D n (4000)) and equivalent width (EW) of the R5 absorption (R5 A ) 
are widely used in AGN's hosts as indicators of the ages of their stellar populations (e.g., 
Heckman et al. 2004; Kauffmann et al. 2003a; Kewley et al. 2006). Incorporating with 
the large high-resolution stellar spectra libraries, the latest theoretical population synthesis 
models of BC03 indicate that both D n (4000) and H<5a are reliable age indicators until a few 
Gyrs after the burst (e.g., Figure 2 in Kauffmann et al. 2003c), although both indices are 
sensitive to metallicity in the very old stellar populations. 

The 4000A break is a strong discontinuity in the optical spectrum of a galaxy. The 
break is mainly caused by the absorption features of ionized metal lines blueward of 4000A. 
The popular definition of D n (4000) is (Balogh et al. 1999; Bruzual 1983) 



r 4ioo f 
J4000 JaQA 

J 3850 jAGtA 



A>(4000) = J , 4 3° 9 ° 5 V;„ (6) 



Gorgas et al. (1999) calibrated D n (4000) in terms of the stellar atmospheric parameters: 
effective temperature, metallicity and surface gravity. In a young stellar population, small 
D n (4000) arises from the low opacity of metal since it is multiply thermal ionized in hot 
stars. On the contrary, large D n (4000) appears in an old stellar population in which the 
metal opacity is expected to be high. 

The strength of the H<5 absorption is an age indicator when the stellar populations are in 
the post-starburst phase, i.e. 0.1-1 Gyrs after the burst when massive stars (O and B stars) 
terminated their evolutions ( Worthey & Ottaviani 1997; Gonzalez Delgado et al. 1999). The 
index H<5a measures the EW of absorption in A type stars, and is defined as (Worthey & 
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Ottaviani 1997) 

E5 A = (4122.25 -4083.50) (1 - (7) 

where Fj is the flux within the AA4083.50 — 4122.25 feature bandpass, and Fc the flux of 
the pseudo-continuum within two defined bandpasses: blue AA4041.60 — 4079.75 and red 
AA4128.50 — 4161.00. As the case of D n (4000), similar calibration has been carried out by 
Worthey & Ottaviani (1997). 

We measure both indices in each starlight spectrum. The starlight spectrum is separated 
from the raw data as described in Section 2.2. For each object, the measured values of 
D n (4000) and H5a is listed in Column (7) and (8) in Table 1, respectively. The comparison 
between the values of _D n (4000) obtained in this paper and that derived by MPA/JHU group 
is shown in the left panel in Figure 2. The same comparison but for H5a is shown in the 
right panel. The left panel shows a significant correlation (with slope ps 1) between the 
two sets of .D ra (4000) values, and the right panel a less significant, but still tight one for 
H^. The comparisons clearly indicate that our measurements are highly consistent with 
those derived by MPA/JHU group, although our B.5a appears to be slightly smaller than 
those of MPA/JHU group at left-bottom corner. This small discrepancy could be ignored 
in current study since we only focus on the statistic properties of the composite AGNs. A 
large discrepancy is only found in one object SDSS J085338. 27+033246.1. The value of H5a 
derived in this paper is 7.22 which is much larger than H5 j4 =2.78 obtained by MPA/JHU 
group. The large discrepancy is due to the difference in the spectra components adopted 
to model both continuum and absorption features. In addition to the seven eigenspectra, a 
powerlaw continuum attributed to AGN is involved in the spectra modeling for the object, 
which reduces its starlight continuum level, and enhances its measured H5a consequently. 



4. Classification of the Composite AGNs 

In this section, we first examine the distributions of the physical properties derived in 
the last section on the traditional BPT diagrams, and then perform a very fine classifica- 
tion on the composite AGNs. Their host's properties are finally explored basing upon the 
classification. 

The BPT diagrams are used as a powerful tool to diagnose dominant energy source 
in emission-line galaxies. Seyferts, LINERs, and star-forming galaxies can be separated 
according to their emission-line ratios. With respect to star-forming galaxies, AGNs have 
larger ratios in [O IIl]/H/3, [N II] /Ha, [S II]/Ha, and [O I] /Ha, because of their stronger 
and harder ionizing field. AGNs therefore concentrate in the up-right corner in each BPT 
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diagram. The theoretical curves discriminating between AGNs and star-forming galaxies in 
the BPT diagrams were recently proposed by Kewley et al. (2001). There are many optical 
classification schemes to segregate LINERs from other emission-line galaxies (e.g., Heckman 
et al. 1980, Veilleux & Osterbrock 1987; Ho et al. 1997). In this study, we adopt the new 
empirical classification scheme recently proposed by Kewley et al. (2006). The classification 
scheme was obtained by analyzing a large sample of emission-line galaxies selected from 
SDSS survey. These new classification line is very similar to that of Heckman et al (1980) 
in both slope and intersection. 

The BPT diagrams are displayed in Figure 3 for the composite AGNs. The three panels 
correspond to the three different methods for classifying galaxies using two pairs of line 
ratios. In each panel, the theoretical classification line is shown by a solid line, and the 
empirical classification line suggested by Kauffmann et al. (2003a) by a dashed line. The 
classification lines used to separate Seyferts from LINERs are shown by dot-dashed lines. 
The size of each point is scaled to be proportional to the corresponding value of _D„(4000) 
simply because of the high level of accuracy of this parameter (See Section 3.2 for details). 

At first, almost all the composite AGNs are located between the theoretical and em- 
pirical lines in the [O III]/H/5 vs. [N II] /Ha diagram, which indicates that our AGN-host 
separations and emission-line measurements are consistent with those of MPA/JHU groups. 
In addition, it is interesting that the objects occupy a narrow vertical stripe with nearly 
constant line ratio [N II]/Ha. The examination of this phenomenon is out of the scope of 
this paper, and will be studied in subsequent works. Moreover, no clear trend is found along 
this stripe for L> n (4000), although the objects with larger D n (4000) are only found at the top 
end of the stripe. In contrast, clear trends are found in the other two panels. The objects 
located in the occupations of LINERs generally have larger D n (4000) with respect to that 
located below the theoretical demarcation lines, which implies that the composite AGNs 
at least could be divided into two sub-groups: H II region-like ionization + young stellar 
population and LINER-like ionization + old stellar population^. 

Since the composite AGNs are initially selected only from the [N II] /Ha vs. [O III]/H/5 
diagram, a more fine classification on the objects is performed as similar as done in Kewley 
et al. (2006). We separate our sample into three groups according to the classification 
scheme described as follows: a) The H II composite AGNs are galaxies that lie below the 
theoretical demarcation curves in all the [N II] /Ha, [S II] /Ha and [O I] /Ha diagrams. These 
galaxies correspond to the "composite galaxies" in Kewley et al. (2006); b) In contrast, 



3 Similar phenomena and indications could be found in the case of L/L^dd and luminosity of [O III 
instead of Z? n (4000), with, however, less significance. 
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the Seyfert+LINERs are galaxies that lie above the theoretical demarcation curves in all the 
three diagrams; c) The remaining are multiply-classified galaxies that are classified differently 
in the three diagrams. These galaxies correspond to the case b "ambiguous galaxies" in 
Kewley et al. (2006). The classification is shown in the three BPT diagrams in Figure 4. 
The demarcation lines are shown in the figure by the same symbols used in Figure 3. Only 74 
objects for which all the line ratios are determined are plotted. The H II composite AGNs, 
Seyfert+LINERs and multiply-classified galaxies are shown by red, blue, and black points, 
respectively. 

Our detailed classification indicates that, among the 74 objects, six ones can be classified 
as Seyfert+LINERs rather than composites. The fraction of these six objects is so small that 
their existence has negligible effect on the overall properties of our sample. Totally, there 
are 45 H II composite AGNs and 23 multiply-classified galaxies. The comparison in the age 
of stellar population between the H II composite AGNs and multiply-classified galaxies is 
illustrated in Figure 5. The dashed line shows distribution of D n (4000) for the H II composite 
AGNs, and the solid line for the multiply-classified galaxies. A narrow distribution with 
median value _D n (4000)=1.30 and standard deviation of 0.11 is found in the H II composite 
AGNs. On the contrary, a very wide and possibly bimodal distribution is identified in the 
multiply-classified galaxies. 



5. Star Formation History and Evolution of AGN 

We plot our sample of the composite AGNs in the -D ra (4000) /H5a plane in Figure 6. The 
dashed line shows the evolution of the standard SSP model with solar metallicity, and the 
dot-dashed line the evolution of the model with continuous star formation history (i.e., an ex- 
ponentially declining SFR ip(t) oc e"*/( 4Gyr )) at solar metallicity (BC03). The D n (4000)/HS A 
plane is a useful diagnosis to determine whether a galaxy has been forming stars continuously 
or in a recent burst. Galaxies with continuous star formation history occupy a narrow locus 
in the plane. A recent burst ended in 0.1-1 Gyrs ago results in a significant displacement 
away from the locus (i.e., be larger in H5a), because the optical spectrum is dominated by 
the emission of A type stars. Even though the small discrepancy in H5a is found between 
our measures and those in MPA/JHU groups, the composite AGNs fall close to the model 
with continuous star formation history, except the object SDSS J085338. 27+033246.1 with 
significant strong H<5 absorption H<5/i=7.22 (see Section 3.2 for explanation). 

Our sample allows us to investigate the issue of co-evolution of AGN and star formation 
in a new approach. Our study differs from the previous investigations in: 1) it should be 
emphasized that each composite AGN certainly contains a Seyfert 1 nucleus in its center, 
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because of the existence of the broad Ha emission. In fact, there might be two populations 
of Seyfert 2 galaxies: with and without a hidden Seyfert 1 nucleus. The spectropolarimetric 
surveys and X-ray observations indicate that the hidden Seyfert 1 nuclei are indeed not 
present in a fraction (10%-50%) of Seyfert 2 galaxies (e.g., Tran 2001, 2003; Moran et al. 
2000, 2001, 2007; Panessa & Bassini 2002; Gallo et al. 2006; Pappa et al. 2001); 2) In 
addition, the basic properties of AGN (i.e., L/L Edd and Mbh) are directly obtained from the 
nuclear emission rather than the tracers with large uncertainties. For instance, the luminosity 
of [O III] emission line is statistically used as a tracer of AGN's bolometric luminosity in 
Seyfert 2 galaxies by many authors (e.g., Kauffmann et al. 2003a; Heckman et al. 2004; 
Kewley et al. 2006). 

It was suspected that L/L Edd might be of significance of "age" of AGN (e.g., Grupe 
et al. 1999; Grupe 2004; Wang et al. 2006). In order to examine this point, the L/L Edd 
vs. -D n (4000) diagram is presented in the Panel A in Figure 7. The diagram indicates that 
the L/LEdd estimated from the broad Ha emission varies as a function of the age of the 
host's stellar population. The black holes with high specific accretion rates tend to reside 
in the hosts with relatively young stellar populations. This dependence therefore strongly 
supports the evolutionary scenario of AGN in which AGNs with high L/L Edd evolve to ones 
with low L/L Edd (Grupe 2004; Wang et al. 2006; Wang & Wei 2006). The Panel C in 
Figure 7 shows a variation of the L/L Edd with the line ratio [O I] /Ha which is sensitive 
to the hardness of the ionizing field. As shown in the diagram, a soft ionizing spectrum 
tends to be associated with a relatively high L/L Edd on average, and a hard spectrum to 
be associated with an extremely low L/L Edd only. This trend is consistent with the strong 
correlation identified between the hardness of hard/soft X-ray emission and L/L Edd for both 
local and distant AGNs (Grupe 2004; Bian 2005; Shemmer et al. 2006). The lowest value 
of the L/L Edd among our sample is 0.003, which is consistent with the prediction of the 
advection dominated accretion flow (ADAF) model. The model can only produce sub- 
Eddington accretions (L/L Edd < 0.01 — 0.1). Because a large fraction of the energy released 
is advected with hot gas into black hole rather than radiated by emission, the gas is extremely 
hot, and produces a harder spectrum in a low state (e.g., Abramowicz et al. 1995; Narayan 
& Yi 1995; Esin et al. 1997; Meyer et al. 2000; Ferreira et al. 2006). 

The another candidate of "age" of AGN is the accretion rate onto a black hole (i.e., 
luminosity released by accretion). Relatively good relations are identified between Ln a and 
-D„(4000), and between L Ha and [O I] /Ha. The two relations are presented in Figure 8. 
Noted that both relations degrade at the right-bottom end. This phenomenon does not 
appear when L/L Edd is used in place of L Ha as shown in the Panel A and C in Figure 7, 
which indicates that the L/L Edd is the major physical driver in the evolution of AGN. 
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A significant correlation between _D n (4000) and the line ratio [O I] /Ha is first shown 
in the Panel B in figure 7. The tight correlation implies that AGN evolves from the high 
state with a soft spectrum to the low state with a hard spectrum. The line ratio [O I] /Ha 
is suggested to be used as a potential indicator of the age of stellar population. 

6. Discussion 

In addition to the AGN's contribution, other main mechanisms of the broad Balmer 
emission lines are the stellar envelopes of Wolf-Rayet stars and OB associations associated 
multiple SN events (Izotov et al. 2007 and references therein). Both mechanisms can produce 
the broad emission lines with gaseous velocity as large as ~ 1000 km s -1 . The luminosity 
of the broad Ha emission line ranges from 10 40 to 10 42 ergs s^ 1 in the composite AGNs, 
which greatly exceeds the typical luminosity of stellar winds from WR or LBV stars (10 36 - 
10 40 ergs s _1 ). Moreover, the broad Ha emission in the composite AGNs is unlikely caused 
by the SN events, because of the existence of the underlying relatively old stellar populations 
in the spectra. 

6.1. Aperture effect 

The redshift ranges roughly from 0.05 to 0.15 in our sample. In this range of distance, 
the typical projected fiber diameter (3") is from ~ 3kpc to ~ 8kpc. When a spectrum 
is taken in the SDSS spectroscopic survey, the fixed fiber aperture would result in more 
contamination of the light emitted from the outer part of host galaxy in distant galaxies 
than in nearby galaxies. First, the effect leads the LINER fraction decreases with redshift 
because the weak emission lines are more difficult to be detected in distant galaxies than in 
nearby galaxies. Such effect exists in our sample, and is shown in the left panel in Figure 
9. The fraction of the LINER-like galaxies with an old stellar population drops significantly 
when redshift z > 0.1, which is consistent with Kauffmann et al. (2003a) and Kewley et 
al. (2006). Secondly, for a fixed aperture, as the contamination increases with redshift, the 
value of -D„(4000) would increase with redshift if a circumnuclear starburst occurs. If so, 
one could see a lower boundary with positive slope in the _D n (4000) vs. z diagram, which 
is, however, not present in the diagram. The aperture effect can not explain why no young 
stellar populations (-D4000 < 1.2 — 1.4) is deteced in the low- redshift galaxies. The negative 
trend shown in the diagram could be explained by the consequence of the increase of L^ a 
with the redshift. Moreover, the lower limit of redshift of z ~ 0.05 ensures that the spectra 
taken within 3" aperture could approximately reflect the global properties (Kewley et al. 
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2005). The right panel in Figure 9 shows the broad-to-narrow Ha line ratio as a function 
of the redshift. As shown by the diagram, in general, our criterion for the existence of the 
broad Ha line is free of the aperture effect. 



6.2. Star formation rate 

In section 4, we have built a sub-sample of the H II composite AGNs which are below 
the theoretical demarcation lines in the three BPT diagrams. The left panel in Figure 10 
shows the [O III]/[0 II] vs. [O I] /Ha diagram for the H II composite AGNs. The distribution 
shows that the H II composite AGNs lie within the occupation of the star-forming galaxies, 
which additionally suggests that hot stars might be main contributor of the ionizing field 
for the narrow emission lines. The luminosity of [O II] emission line is an indicator of 
current star formation rate (SFR) for H II regions as suggested by Kennicutt (1998) and 
Kewley et al. (2004). However, the origin of the [O II] emission is complex in the H II 
composite AGNs, since both AGN and starburst can contribute to [O II] (e.g., Kim et al. 
2006; Yan et al. 2006). In order to estimate the contribution by AGN, the right panel 
in Figure 10 shows the luminosity of [O III] as a function of the line ration [O II] /[O III]. 
The solid line is the least-squares bisector regression for Type I AGNs given by Kim et al. 
(2006). The H II composite AGNs clearly show enhanced [O II] /[O III] ratios. Similar as 
done in Kim et al. (2006), we estimate the [O II] luminosities emitted from H II region by 
assuming the enhanced [O II] /[O III] ratios are caused by star formations. The mean value 
of the [O II] /[O III] ratio is ~ 1.0 for the H II composite AGNs. Given the average [O III] 
luminosity L[om] ~ 3 x 10 41 erg s -1 , the [O II] /[O III] ratio predicted by the regression is 
only ~ 0.25. Then, on average, ~75 percent of the [O II] emission is estimated to come from 
star formation in the H II composite AGNs. 

We use the recent calibration (Kewley et al. 2004) 

SFR - 7 - 9X 16.73 - 1 7^"grO/H) + 1 21 M ^' 1 < 8 > 

to estimate current SFRs, where L[oii],42 is the luminosity of [Oil] emission in units of 
10 42 ergs s -1 . The metallicity is fixed to be log(0/H) + 12 = 9.2, i.e., twice of the solar, as 
suggested by Ho (2005b). 

A correlation between the current SFRs and L/L E dd is shown in the left panel in Figure 
11. The relatively large scatter could be caused by both the underestimation of L^ a and the 
contribution of AGN to L^qu]- The correlation indicates that a more powerful AGN tends to 
be associated with a more intense current star formation. The relation is reasonable because 



-17- 



the strengths of both AGN and star-forming activities depend on the consumption of gas. 
Heckman et al. (2004) found a roughly constant ratio of the star formation rate to the accre- 
tion rate onto black hole in the SDSS survey. The dependence of the SFR on the luminosity 
of AGN is further supported by the correlation, between the X-ray luminosity in 0.5-4.5 keV 
band and CO luminosity, found in Seyfert 1 galaxies and in quasars (Yamada 1994). The 
middle and right panels in Figure 11 shows variations of the SFRs with D n (4000) and B.5a, 
respectively. As shown by the diagrams, the current SFR decreases when stellar populations 
passively evolves, because of the depletion of the gas. Both relations are consistent with 
Figure 6 in which a continuous star formation history with an exponentially decaying SFR 
is required to match the observations to the modeled evolution locus. In addition, the [O II] 
inferred SFRs are typically below lOM yr _1 . The suppressed current SFRs are identified in 
a number of luminous AGNs by several authors (e.g., Ho 2005b; Kim et al. 2006; Wang & 
Wei 2006; Martin et al. 2007). 



6.3. Role of L/L E dd in AGN evolution 

Similar as the scenario proposed by Weedman (1983), Wang & Wei (2006) recently 
proposed a scenario of co-evolution of AGN and star formation in its host basing upon 
the differential growth of the black hole mass and bulge mass. In the scenario, after a 
starburst occurs first and triggers accretion onto the central black hole latter, AGN then 
shines with decreasing L/L^dd and insignificant increasing black hole mass, as the young 
stellar population ages and fades. The scenario is supported by the prediction of the model 
developed by Kawakatu et al. (2003). The model predicts that the AGN-dominated phase 
follows the host-dominated phase in a time scale of a few 10 8 years. This change in phase is 
also predicted in the more solid physical models developed by Granato et al. (2001, 2004). 

In addition to the models, there were accumulating observations supporting that L/Lsdd 
plays an important role in the evolution of galaxy in recent years. Mathur (2000) suggested 
that NLS1 might be young AGN with high L/L Edd and small black hole mass. The young 
stellar populations within post-starburst phase are found to be associated with NLSls in a 
few cases (e.g., Wang et al. 2004; Wang & Wei 2006; Zhou et al. 2005). The evolutionary 
significance of L/L^dd was suspected by Grupe (2004). In two years latter, Wang et al. 
(2006) extended the well-documented El space into the infrared color a(60, 25) that is an 
indicator of relative strength of star formation with respect to AGN. The infrared-color 
dominated El space implies a relation between the L/L^dd and star formation history. With 
the advent of SDSS, Kewley et al. (2006) and Wild et al. (2007, and references therein) 
revealed a relationship between the star formation history and I/[oni]/c* in type II AGNs . 
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Taking advantage of the composite AGNs, the evolutionary significance of L/L E dd (or 
the El space) is obviously demonstrated by the clear dependence of L/L Edd on _D„(4000). In 
addition, the evolutionary significance of the El space is also implied by the tight correlation 
between [O I] /Ha and _D n (4000). The line ratio [O I] /Ha is known to be sensitive to the 
hardness of the ionizing radiation field. The X-ray photon index is an important component 
in the El space (e.g., Wang et al. 1996; Sulentic et al. 2000; Xu et al. 2003; Grupe et 
al. 1999; Grupe 2004). In summary, by adopting the new approach, our results are not 
only consistent with the previous studies, but also strongly support the argument that AGN 
evolves along the El sequence from high L/L-^dd end to low L/L-^dd end as stars evolve from 
young to old stellar population. 



6.4. "Young" and "Old" composite AGNs 

Our results indicate that ~80% of the composite AGNs are high L/L Edd emitters asso- 
ciated with a post-starburst (with age less than 1 Gyr, corresponding to a value _D n (4000) < 
1.5). In the context of the evolutionary scenario, the composite AGNs associated with both 
high L/L E dd and young stellar population might be "Young" AGNs. The spatial connection 
between the post-starburst and AGN is supported by the spatial resolved spectroscopy in 
the three AGNs selected from SDSS (Goto 2006), and by recent photometric observations 
with spatial resolution of 0.1-0.2" (Davies et al. 2007). In terms of the equivalent width of 
the narrow Ha emission, Westoby et al. (2007) found the composite galaxies peak in the 
valley between the red and blue sequence (e.g., Strateva et al. 2001; Baldry et al. 2006), 
which suggests the composite galaxies are in the transition state. By calculating a series of 
photoionization models, Stasinska et al. (2006) suggested that the right wing in the BPT 
diagram could be explained by the variation of the balance between massive stars and AGN 
powers along the wing. In the BPT diagram, a galaxy moves toward right-up direction from 
the occupation of the composite AGNs as the fraction of AGN's power increases. 

On the other hand, although the multiply-classified galaxies mostly lie below the theo- 
retical demarcation curve in the [N II] /Ha vs. [O IIl]/H/3 diagram, a fraction of them are 
located in the LINER region in the [O I] /Ha vs. [O III]/H/3 diagram. The harder ionizing 
field and smaller AGN's power could be inferred in these multiply-classified galaxies, because 
of the saturation of [N II] /Ha with increasing nebular metallicity in H II region (e.g., Kew- 
ley & Dopita 2002; Pettini & Pagel 2004). We suspect that the multiply-classified galaxies 
associated with very old stellar populations are at the end of life of AGNs. 



-19- 



7. Summary 

The goal of this paper is to study the evolutionary connection between AGN and its host 
galaxy in partially obscured AGNs (i.e., Seyfert 1.8 and 1.9 galaxies). Not as the previous 
studies, we can derive the properties of both AGN and its host together in these objects by 
direct measurements rather than some tracers. Taking advantage of the SDSS spectroscopic 
survey with large number of spectra of emission-line galaxies, we searched for Seyfert 1.8 
and 1.9 galaxies (composite AGNs) from the emission-line galaxies (MPA/JHU catalogs), 
which are located between the empirical and theoretical classification lines in the [O IIIj/H/3 
vs. [N II] /Ha diagram. The analysis of these 85 composite AGNs allows us to draw results 
as follows: 

1. The distribution of the composite AGNs on the [N II] /Ha vs. [O IIl]/H/3 diagram 
shows a narrow vertical stripe with nearly constant line ratio [N II] /Ha. In contrast, 
the composite AGNs can be separated in the other two BPT diagrams, in which the 
objects associated with the very old stellar populations are located in the occupation 
of LINERs. 

2. The BPT diagrams allow us to further divide the 68 composite AGNs into two sub- 
groups: 45 HII composite AGNs which are below the theoretical separation lines in 
all the three BPT diagrams; 23 multiply-classified galaxies which are classified differ- 
ently in the three BPT diagrams. All the HII composite AGNs are associated with the 
young stellar populations (<lGyrs), while a very wide and possibly bimodal distribu- 
tion of the age of the stellar population is identified in the multiply-classified galaxies. 
Totally, ~80% of the composite AGNs are associated the young stellar populations 
within post-starburst phase (L> n (4000) < 1.5). We suspect that the composite AGNs 
associated with the very old stellar populations are at the end of life of AGNs. 

3. The broad Ha emission inferred L/LEdd forms a smooth sequence with the age of the 
stellar population of AGN's host as assessed by L> n (4000), which confirms the previous 
studies using other tracers (e.g., £[om]/c*, infrared color a(60,25)). The sequence 
degrades when L/L^dd is replaced by Lh«. In addition, the L/LEdd roughly depends 
on the line ratio [O I] /Ha which is sensitive to the hardness of the ionizing field. Both 
theoretical population synthesis models and [O II] inferred SFR ~ 1 — 10M Q yr _1 show 
a suppressed current star formation in the composite AGNs. These results strongly 
support the evolutionary scenario suggested in Wang & Wei (2006). 

4. We first find a significant correlation between the line ratio [O I] /Ha and L>„(4000) in 
the composite AGNs. If it is confirmed, the line ratio could be used to trace the age 
of stellar population in type I AGNs. 
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Fig. 1. — Left column: the modelings and subtractions of the starlight components by the 
linear combination of the seven reddened eigenspectra in the three typical cases. In each 
panel, we plot the emission-line spectrum, observed spectrum and modeled starlight spectrum 
from top to bottom. The spectra are vertically shifted by arbitrary amounts for visibility; 
Middle column: the modelings of the emission-line profile by a set of Gaussian components 
for H/3 region; Right column: the same as the middle column but for Ha region, [see the 
electronic edition of the Journal for a color version of this figure.} 
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Fig. 2. — The comparisons between the value obtained in this paper and those derived by 
MPA/JHU groups for both indices -D n (4000) (left panel) and H5a (right panel). It is clear 
that the two independent measurements are highly consistent with each other. 
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Fig. 3. — The three diagnostic BPT diagrams for the composite AGNs. The theoretic 
demarcation lines separating AGNs from star-forming galaxies proposed by Kewley et al. 
(2001) are shown by the solid lines, and the empirical line proposed by Kauffmann et al. 
(2003a) by the dashed line. The dot-dashed lines drawn in the [S II] /Ha vs. [O IIl]/H/3 and 
[O I] /Ha vs. [O IIl]/H/3 diagrams show the empirical separation scheme of LINERs proposed 
in Kewley et al. (2006). The size of each point is proportional to the corresponding value 
of £> n (4000). In the [N II]/Ha vs. [O IIl]/H/3 diagram, the composite AGNs distribute in 
a narrow vertical stripe with nearly constant line ratio [N II] /Ha. However, the composite 
AGNs with very old stellar populations (large D n (40000) are mostly located in the LINER 
region in both [S II] /Ha vs. [O IIl]/H/3 and [O I] /Ha vs. [O IIl]/H/3 diagrams. 



-29- 




Fig. 4. — The three diagnostic BPT diagrams showing our fine classification on the 74 
composite AGNs. The displayed curves are the same as in Figure 3. The HII composite 
AGNs located below the theoretic lines in all the three BPT diagrams are presented by the 
red solid squares, and the Seyferts+ LINERs located above the lines by the blue circles. The 
black open squares represent the multiply-classified galaxies which are classified differently 
in the three diagrams, [see the electronic edition of the Journal for a color version of this 
figure] 
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Fig. 5. — The histograms of the index D ra (4000) for the HII composite AGNs (the dashed 
line) and for the multiply-classified galaxies (the solid line). All the HII composite AGNs 
are associated with young stellar populations (_D n (4000) < 1.6), while the multiply-classified 
galaxies show very broad and possibly bimodal distribution. 
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Fig. 6. — The D„(4000)/H5a plane for the composite AGNs. The dashed line represents the 
stellar population evolution locus of the SSP model with solar metallicity, and the dot-dashed 
line the model with exponentially decreasing star formation rate ip(t) oc e~'^ 4Gy \ The model 
with continuous star formation history is required to explain the composite AGNs. 
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Fig. 7. — Panel A: A sequence of L/L-^dd with the age of stellar population of host for the 
composite AGNs. Panel B: A tight correlation between D„(4000) and the line ratio [O I] /Ha. 
Panel C: L/L-^dd plotted against the ratio [O I] /Ha. As shown by the plot, there seems to 
be a trend that larger [O I] /Ha is statistically associated with lower L/L-^dd- Panel D: A 
tight correlation between the luminosity of the Ha broad component and that of [O III]. The 
best unweighted fitting is represented by the solid line. 
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Fig. 8. — The luminosity of the Ha broad component plotted against -D n (4000) (left panel) 
and RSa (right panel). Although L Hq is related with both D„(4000) and US a at the high 
luminosity end, both correlations degrade at the low luminosity end. 
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Fig. 9. — The parameter _D„(4000) (left panel) and broad-to-narrow Ha line ratio (right 
panel) as a function of redshift. 



-35 - 




Fig. 10. — Left panel: The [O III]/[0 II] vs. [O I] /Ha diagnostic diagram for the composite 
AGNs. The solid lines represent the classification scheme proposed by Kewley et al. (2006), 
and the dashed line the LINER separation line in Heckman (1980). Right panel: The line 
ratio [O II] /[O III] plotted against the luminosity of [O III]. The solid line shows the least- 
square regression for type I AGNs given by Kim et al. (2006). The composite AGNs clearly 
show enhanced line ratio [O II] /[O III]. 
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Fig. 11. — [The O II] luminosity inferred SFR plotted against LjL^A {left panel), _D ra (4000) 
(middle panel), and H5a (right panel) for the composite AGNs. 
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Table 1 — Continued 
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i^n(,4UUU) 


MA , 
ilOA 

( A ) 


\ W) erg s ) 


lOg(M / Mq ) 




t i\ci TTT^ 
L([U ill]) 

erg s ) 


(>-) 


( 2 ) 


V>) 




(5) 


w 


(<) 


(°) 




(,10) 




(12) 


SDSS J151732. 81+602844. 6 


0.129 


0.86 


0.37 


0.075 


7.71 


1.47 


-0.240 


0.68 


7.17 


0.08 


45.87 


SDSS J152023.49+502947.8 


0.093 


0.42 


0.31 


0.026 


0.59 


1.28 


3.113 


0.43 


7.58 


0.02 


7.48 


SDSS J152036. 05+312225. 8 


0.107 


0.51 


0.31 


0.047 


2.29 


1.44 


0.943 


0.94 


8.05 


0.01 


21.32 


SDSS J153015.55+531035.1 


0.174 


0.43 


0.26 


0.035 


1.85 


1.24 


4.034 


1.39 


7.56 


0.06 


62.76 


SDSS J163330. 21+213634. 9 


0.095 


0.52 


0.57 


0.176 


2.28 


1.84 


-2.150 


0.18 


7.28 


0.02 


12.14 


SDSS J164643.66+201654.5 


0.081 


0.64 


0.36 


0.048 


1.11 


1.46 


1.113 


0.27 


7.34 


0.02 


2.65 


SDSS J163413. 74+445947. 8 


0.058 


0.55 


0.38 


0.058 


2.46 


1.38 


1.427 


0.17 


6.83 


0.05 


8.65 


SDSS J164531. 50+403224.6 


0.112 


0.50 


0.31 


0.051 


1.99 


1.26 


3.600 


0.32 


7.16 


0.04 


32.95 


SDSS J164552.20+442614.8 


0.070 


1.37 


1.06 


0.511 


5.74 


1.50 


-1.040 


0.61 






5.04 


SDSS J165039.15+634429.0 


0.148 


0.42 


0.27 


0.036 


1.42 


1.22 


3.390 


0.50 


7.57 


0.02 


25.31 


SDSS J171855.43+275947.4 


0.111 


0.49 


0.28 


0.039 


1.34 


1.40 


2.684 


0.36 


7.13 


0.05 


11.94 


SDSS J210641. 11-075501. 2 


0.131 


0.64 


0.48 


0.080 


1.13 


1.38 


2.289 


0.15 


6.29 


0.17 


9.77 


SDSS J210721. 91+110359.0 


0.042 


0.26 


0.47 


0.103 


2.52 


1.47 


0.996 


0.06 


6.55 


0.04 


W 1.31 
0.72 


SDSS J213741.44+005316.9 


0.041 


1.16 




0.081 


3.23 


1.58 


-1.180 


0.04 


6.50 


0.03 


SDSS J220531. 56-073734.8 


0.060 


0.85 


0.77 


0.161 


3.18 


1.41 


-0.500 


0.09 


7.35 


0.01 


0.73 


SDSS J222705.69-090424.8 


0.073 


0.49 


0.30 


0.032 


1.25 


1.45 


0.568 


0.10 


6.63 


0.05 


1.36 


SDSS J225057.27-085410.9 


0.065 


0.45 


0.41 


0.045 


2.54 


1.28 


2.756 


0.08 


7.00 


0.02 


17.22 


SDSS J225102.65-093824.1 


0.081 


0.40 


0.40 


0.036 


0.65 


1.37 


2.976 


0.09 


6.81 


0.03 


2.42 


SDSS J232643.08-095927.8 


0.070 


0.42 


0.36 


0.052 


4.03 


1.25 


-0.22 


0.13 


6.51 


0.09 


6.00 


SDSS J232714.52-102317.5 


0.065 


0.52 


0.30 


0.028 


0.95 


1.23 


0.991 


0.18 


6.55 


0.11 


7.92 


SDSS J233942.54+142305.3 


0.149 


0.49 


0.22 


0.033 


1.69 


1.21 


2.105 


1.69 


6.95 


0.29 


79.06 



a The spectra ID is spSpec5 1608-0267-405 given by SDSS. 
b Thc spectra ID is spSpec51989-0480-120 given by SDSS. 



